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Abstract 

In this paper, the photon stationary transport equation has been extended from R 3 to C 3 . A solution 
of the inverse problem is obtained on a hyper-sphere and a hyper-cylinder as X-ray and Radon transform, 
respectively. We show that these results can be transformed into each other and they agree with known 
results. 
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Introduction 


The transport equation has many applications in different fields of science and technology 
(see 1]- 7] and the references therein). The solution of the stationary non-diffusive transport 
equation for photons solves many imaging inverse problems such as Computerized Tomog¬ 
raphy (CT), Single Photon Emission Computerized Tomography (SPECT) and Positron 


Emission Tomography (PET) a -llQl 


In three dimensions, the solution of transport equation, without attenuation effect, ap¬ 
pears as the “X-ray cone beam” transform, without any restriction on the set of source points 


x m 


11 


-13, [15]. The reconstruction formula is obtained by averaging the X-ray data on the 
unit sphere in R 3 . A special case with application to tomography considers point sources 


lying on a space curve and is treated in 


U 


16-LL8J. B D Smith has introduced an approach 


for converting divergent beam data into parallel beam data in order to use the well-known 


Radon inversion procedure [19j. In another work, P Grangeat has made a conversion of 
X-ray data into 3-dimensional Radon data in order to apply the Radon inverse formula 20]. 

In two dimensions, Novikov jg] has worked out an explicit inversion formula for transport 
equation where attenuation effects are taken into account. In his work, Novikov has used 
an extension R 2 to C and apply the procedure of analytical continuation. The inversion 
iroblem of the transport equation with attenuation in three dimensions has been solved by 


m using quaternion analysis. The aim of this paper is to use several-complex analysis to 
find an inversion formula for the transport equation in R 3 . We would also like to use this 
idea to find an explicit relation between the x-ray transform and Radon transform solutions. 
We will also compare our results with existing ones in the literature and we show that they 
are the same. 

This paper organized as follows: in the next section we complexify the photon transport 
equation. Then, in section 2 we obtain the inverse problem solution for the unit ball subspace 
as x-ray divergent beam transform. In section 3 we use cylindrical coordinate subspace to 
find the solution. After that we present the results for the three-dimensional case in section 
4. In section 5, we will convert the two results into each other. Finally a conclusion and 
perspectives are given in the last section. 


2 


























I. PHOTON TRANSPORT EQUATION IN SEVERAL COMPLEX 

Photon transport is considered as non-interacting phenomenon with light speed, c. For 
simplicity we set c = 1 and consider the stationary photon transport equation 


n ■ V x $i(x) = pi (x), 


( 1 ) 


where $i(x) is the photon flux density, pi £ Cq the photon source term and n = (n i, n 2 , n 3 ) 
is the photon motion direction. We now consider the analytic continuation of this equation 
in several complex variable space. In two dimensions, this procedure is known, i.e. x —)■ 
x + iy , (y —>■ 0). In three dimensions, analytic continuation requires complexification of 
the photon direction vector, as done in !§, 9, 21], This can be done by adding a second 
photon transport equation with respect to a variable y in a perpendicular direction to the 
initial photon direction, n. We call this direction = (n^,n^,n^). For spatial dimension 
equal or greater than 3, the definition of rr 1 is not unique: can be any vector in the 

perpendicular hyper-plane to n . We do not worry about definition of because we will 
see that in the final step we take an average over all possible directions. Thus, this procedure 
leads to the following equation 


(n-V a , + n ± -V J/ )$(x,y),= p(x,y) (2) 

where p(x,y) is the source term extending pi(x). The above partial differential equation is 
defined in an R 6 = (R 3 x R 3 ) Euclidean space of (x, y). We now embed it into C 3 . We first 
define 

' 

r = x + i y 

< V = V x -iV„ (3) 

N = n + i n x 

and then we impose the following condition 

( n_L • V x — n • Vy) <f>(x, y) = 0, (4) 

where $ is a function defined on a domain in a closed set Q C R 6 , this condition means 
that $1^ = 0, where S = {zi £ C 3 ,|zj| = 1, i = 1,2,3} is the Shilov boundary. The above 
condition means p is a function with compact support on dtt. This condition also means that 
there is a holomorphic function such as 0 £ Q so that dcj) = 0 and (f) has an holomorphical 
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extension on 22], such as $. Now, we can write the equivalent of ([2]) in C 3 as: 


N • V<h(x, y) = p(x, y), (x, y) G ft C 

$ls = 0, 


( 5 ) 


In fact, the above equation is a simplified stationary transport equation with source or sink 


term and without attenuation 


divergent beam integration 32] 


23]. Since N ■ V is a directional derivative, its inverse is x-ray 


(fbp)(r,N)= / dt p {r + tN). 


( 6 ) 


The physical boundary condition on (fhp)(r,N) is the following. For a given direction 
N, because of the support hypothesis and the prescription on the direction of integration, 
(dfp)(r, N)|g = 0, and N points outward of dXl. 

Here N depends only on three independent parameters because of n • rr 1 = 0, = 


In 3 -1 = |n| = 1 and also iN = N 3 - which yields: N • N 3- = 2i. 


The x-ray divergence beam transform of p is a superposition of Radon transform in 

and in (—N)-direction ((7£p)(r, — N) = 
27]: 


N-direction for ((Jlp)(r,N) = f R dt p(r + fN 
Ldt p(r-fN)), i.e. x-ray transform of p is [2lJ, 


(*p)(r,N) = (Kp){ r,N) + (Kp){ r,-N) 


(7) 


In fact, the above relation and the Paley-Wienner’s theorem for the Radon transform 24] 
allowed us to extend compactly supported space to a subspace of C 3 with an holomorphical 
continuation. For fully complexified equation ([5]) We have to extend N in complex form also. 


This extension will be done by generalizing the method of 
parameters and variables: 


3 , 21 ]. 


Let us define the following 


A* = | Xi | e l9i = \Xi\(rii + in 3 -), 

Zi = Xi + hn, 

zi = Xi — i i/i, i = 1,2, 3. 


( 8 ) 


where is the angle between o^-axis and x following Lemma gives the holomorphical ex¬ 
tension of equation ([2]). 
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Lemma 1. The extended several complex variable form of equation (U|) is: 

3 


VyKr,A).£^M = p(r) 


( 9 ) 


2=1 


where rj is defined in m 

Proof. By choosing the parameters and variables in the equation (|8j) . the derivative 
operator in ([2]) becomes as 

3 


_ , _ x 1 \ - w d Id 

n.V I + n • S - + - 57S - 

i =1 x 


( 10 ) 


Extending the above equation to complex, A*, yields 


N 


N 


(n ■ V* + n 2 -• V t )-> V* = = - £ |A 


2=1 


2=1 

N 


JOi 


d 


+ 


d 


dzi | Ai| e Wi dli 


1 f . d 1 d \ 

2 ^ V l ® z i + Kdzi' 
2=1 x 


( 11 ) 


Here Va is the several complex form of left-hand side of (ITU and when |Aj| —> 1 then, Va —> 
‘left-hand side of (TTO]) ’. Now consider the following change of variables 


(/ i . (Aj Zi Aj zfj . 


Setting r/i = 3FJ(?7i) + i^s(rji), we have 


£(7/0 = (|Ai| + lAil” 1 )^, % Vi ) = (|Ai| - |Ai| _1 )2i, 


( 12 ) 


(13) 


where 


(14) 


U = (r • N)j = Xi cos 9i + yi sin 6 U 
Si = (r • N^); = -Xi sin 0* + yi cos 0* , (i = 1, 2, 3) 
where r = (x, y) = (xi,x 2 ,x 3 ,y 1 ,y 2 ,y 3 ) is a vector in M 6 , ( 6 C 3 , t e (U,^,^) and 


s = (si,s 2 ,s 3 ). We can write r = X)(=i (^ L + Nf~ sf). We also note that x = n;Xj 

and y = nf yi, so that we have x ■ rU = 0 and y • n = 0. Using the change of variables 

introduced in (TT2lh the Va* of equation (fill) becomes 


v ^K A 4 + yl;)^U' 2 -w 2 ) w 


Dehning now the function 


<Wv) = MSKriMv)) = (|A,| 2 - IA,| -2 ) 4>i(r, A), 


(15) 


(16) 
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we see that 


(17) 


Va 8 *i(r,A) = 


d$j(i 7 ) 
drjl 


Hence the several complex variable extended form of (J2J) is obtained as equation 


II. HYPER-SPHERICAL COORDINATES 


To solve equation ()9|) . we shall use the following theorem in the several complex analysis. 


Theorem 2 


26 


28]. If h(z) = Y^i=i^ l i{ r i)d , q i is a compactly supported (0,l)-form in C N 


of class C 1 that is d-closed. Then, the solution of this problem is known as: 


f{z) = f h(Z)AB(£,z) (18) 

Jc N 

then f is a compactly supported function such that df = h. Here B(£,z) is the Bochner- 
Martinelli kernel which is defined as follows: 


B&z) 


(N — 1)! 
(27ri) w 


N 




(£i z i) 


N 

dw i (Z) = ' 52 (-l) i+ 1 G i (Z,z) 

2=1 




(19) 


where 


dwftf) = d£ 1 A • • • A dfi-i A dfi + 1 A • • • A d£jv A d£i • • • A d^.A (20) 


We rewrite now @7]) in cM> c = ^ dip form as follows 

d$c(v) = hfrf) (21) 

where h(rj) = Piiv) dr/, : • Thus, in equation fl2T]) h is a one-form function with compact 

support in il C C 3 . Pi(v) is the several complex variable form of p and in our case Pi{rj) is 
the same for all i. We should add that p (equivalently pfi])) is a holomorphic function, i.e. 

= 0, Vz G {i\i = 1, 2, 3}, thus, it satisfies the integrability condition in several complex 
variable analysis fp- = Now, we use the above theorem to obtain the solution of our 
system of equations. But we must prove first the following Lemma. 

Lemma 3. The solution of the inverse problem equation m is: 

^)=^£/^4 ( * p)(r ’ N)dr (22) 

2=1 J 1 
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with Ni-§= = (m + inf) [£: + i^-). 

Proof. By using Theorem 2 we can write <f>i as: 

iiM= f (|A,| 2 -|A ( |- 2 ) _I ft® GiK -r,)ds 
J n 

where df indicates the integration volume element defined as 

di = (-l) i+1 dli A dwi{£) = d& A d£i • • ■ A d£ 3 . 


Then 4>i in equation 


can be written as: 


®i(v) = 


6 - Vi 


PiiO 


( 23 ) 


(24) 


(25) 


4tt 3 Jn If - v\ 6 

Setting |A| —* 0 and making the following change of variables (£,- — rji) —> r/ ? ; = A' 1 z* = 
A£ x (r • N)j = A” 1 Ni t means that we have considered a source in £ with photon flying in 
the N—direction. This means that r* —* r* + A£ A Thus, by using equation (Tlfil) the above 
formula yields 4>(r, N) as follows 


tva 1 f A/,; t p(r + N f) / |. 12 1 . | 

®( r " N l = 7^3 _/ n 2 3|A|-6i6 |A,|-2 |W=EN 


(26) 


2=1 


Using | Aj| A£ = A* and spherical volume element for d£ = 2 3 |A| 6 f 5 dt dT, t e M + , we obtain 


3 

4(r,N)= AV f A, A) f'.Vpn'r. N) 

ln 


dT 


(27) 


For | Ai| = 0, we have also Va = | \ an d equation (|22|) is proved. A 


III. HYPER-CYLINDRICAL COORDINATES 


To obtain the solution of the equation Q in the cylindrical coordinates we shall use the 
following theorem: 


Theorem 4 29]. The fundamental solution of the following problem in C N 

d 


dz 


G(z) = <5(x,y), z(= x + iy) G C ,x,yG 


t N 


is: 


N 


G(z) = £ G,(z) = - 


2=1 


(TV - 1)? 

(2iri) w 


N 

E 

2=1 


127V ' 


(28) 


(29) 
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We prove first the lemma. 


Lemma 5. The solution of the equation 

VaG(ij) = <5(t,s), (i = 1,2,3) 


(30) 


is: 


GM = E5m. <AW = - sgn(|Ai l 2 ~ |Ai| 2) 

z ' Z 7r° 

i=l 


23 Vi 


M r 


(31) 


Proof. We start with X7\G(v) = A(9ft(?7), ^(77)), ?? 6 C 3 and use equation (TT5T) . Hence 


1 

2i 


I 2 1) 

,1-2) 9 | 

f - 1 

1 r 


^4^1 \v \ 6 J 


= 6{3t{rH)Mrk)) 

= 2 -2 | |Aj| 2 — |Aj| -2 | 5(s, t), 


(32) 


which, after some algebraic steps, leads to the equation (j3T]) . A 

For |Aj| = (1 ± e), e —> 0 + , we have 
T1 A(^ 


GfM = 


2 7r 3 2 5 s 4 (s,- ± i ets 


2 6 7T 3 


= = Ft^—T ( P -V.^— =F ivrsgn(ti) 5{si) 

S Si 


(33) 


Now, the Green function introduced in the above equation allows us to construct the limiting 
value <f> of the solution of 

Va 4>(r, A) = p(r) (34) 


Theorem 6. The solution of the inverse problem for the equation |) is obtained as 


"W=-^i:/^ 1 ^(Wr,N)< i r. 

i—1 


(35) 


Proof. The solution <!>, is given by 


^(r,A)= / Gi(r — r', N) p(r') dV' — / Gtf, N) p(r - r') dV'. 
Jn J o 


(36) 


The two solutions for A —)■ 1 ± 0 + becomes <f> ± 


as: 


*f(r,N) = / Gt(vfN)p(r- r ')dV' 
Jn 


(37) 
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Substituting Gf and dV' = 2 6 dt' s' A ds'dT 5 as hyper-cylindrical coordinates (dr 5 is the 
spherical surface element of 5-sphere) we obtain 

$f(r,N)= / (n i n P ){r,N)dT 5 + ^ [ [ p{r - Nt') sgn(f') dt' dV 5 . (38) 

^ J r 5 tt O J R 


Using 0 = <f> + — <f> we obtain 


0i(r,N) =-- {(KiTlp)( r,N)) 5 


7T 


(39) 


where ([•])*, = f Fk [•] dr*. indicates average on the k-sphere (in our notation we eliminate 
subscript for k = 6 such as V = T 6 ), and "H* is the Hilbert transform over Sj and is defined 
as 

1 r f(r i 

(40) 


Hif : —P V. [ ^dX; 

7T ./ 


where ‘P.V.’ indicates the principal value of a singular integral. Finally, this result obtained 
in cylindric coordinates should be consider as in a cylinder of C 3 . Using Cauchy integral on 
the variable, which is not on the unit ball (to project on Shilov boundary) Lemma on page 


32 of 


22] gives 


A±L, 


,1=1 di Aj 


dfa, (A G o \ r) 


(41) 


where <f> is a holomorphic extension of 0 = <f> + —<f> . For A, —y 0 we have Va = | ]y0=i \ 1 ~§Y 
and 0. = — (l + — + •••). Thus, the following formulas hold 

x Mi \ Mi J 

(42) 

(43) 


' w - -AsAL, 


<f>i{ t,s) 

2 


2—1 

3 


— Y — [ 

1 1 U J\H 1 = 1 


1/4* l = i /L 

/r- 1 (H^p(t,s)) 5 ^, (A e n\r) 


On the unit ball = F, it gives equation (j55j) . ▲ 


IV. THREE-DIMENSIONAL RESULTS 


When y —> 0, g(r) —> g(x) then V y g(x) = 0 and dr —A 4-71 dr 3 where T 3 is the surface of 


a 3-sphere and by using (J4]), we also have , N t 


-l a_ 
* cfTi 


n ■ V x , Thus, 


p(x) = 


2 7T 3 


n • V, 


'r 3 


-2tt Ell n) 

(Tp)(x, n) 


dr. 


x e 


(44) 
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This equation includes the x-ray transform, a result obtained by 10]. Moreover we can also 
show that the term containing the Radon transform is equal to the other term. 

Remark 7. g = ('H i Tlp) 5 is an average on a 5—sphere, while g — (X p) is an average over 
a 6—sphere. In reality, this lack of averaging over an angle in the cylinder is compensated 
by an holomorphical extension via Cauchy integral. In analogy with potential theory, the 
Cauchy kernel is the same as ?’ -1 potential, which is a central field potential in physics and it 
describes circular motion of a particle. On the other hand, the Cauchy integral in equation 
(|4T| averages over free angles and completes the averaging over a 6—sphere. This means 
that the 6—cylinder is topologically converted to a 6—sphere. In tomography, this procedure 
fulfills Tuy’s x-ray inversion formula condition, which states the “For any bounded object, a 
curve consisting of two circles [...] or some type of spiral curve will satisfy all the conditions 


on the curve required in the theorem.’ 


14] In three dimensions, Tuy’s condition is fulfilled 


by a sphere because great circles exist on a sphere. 


V. CONNECTION BETWEEN TWO APPROACHES 


We also verify the inter-convertibility these results by using the following proposition: 


Proposition 8 30]. Suppose that f G Cq and x6i 3 and 0 G § 2 . Then 

[ Xf(x, • ni)dT 3 = -(TZf)'(x • n).A 

J r 3 


(45) 


Lemma 9. Using the above proposition we have (Xp) = —2 tt {TLiJZp). 
Proof. Here some work on the left-hand side of this proposition yields 


Xf(x,n 1 )5\n-n 1 )dT 3 = - (Xf)'(x, ni) d(n • ni) dT 3 


'r 3 


'r 3 


f (X J)'(x. ni) 
r 3 (n • ni) + iO 
2 7r TL(X /)'(x, n), 


+ 


'r 3 


W/(x, 

(n ■ ni) — iO 


dr, 


(46) 

(47) 


where ‘(•) / ’ is the differentiation with respect to m. In equation (1461) . we have used the 
definition of 5 —function as a generalized function [311]. Note that in 'H(Xf)' the Hilbert 
transform is relative to n. Equating this result to the right-hand side of the equation in the 


10 








previous proposition, we obtain 2 n) = (7?./) , (x- n). By taking another Hilbert 
transform on both sides of this equation and using “H 1 2 = — V we obtain —2 7r (Xf)'(x, n) = 
n(nf)\x- n). These two last relations between the x-ray divergent beam transform and the 
Radon transform, equation (13) imply that — 2 it (A/)(x, n) = ( "H7^/)(x • n). The following 
equation holds: ( 7/7?./) (x ■ n) = J^ =1 ('H.,:7^/)(x, n), where Hi is Hilbert transform relative 
to Si . This relation finalizes our proof. ▲ 


Conclusion 

A solution of the transport equation in three dimensions without attenuation effect is 
obtained by using its extension to the several complex variable space C 3 4 5 6 7 8 9 . We have shown that 
the solution of this extended transport equation in two different subspaces with cylindrical 
and with spherical coordinates could give rise to two different approaches of solution in terms 
of Radon and x-ray divergent beam transforms respectively. However, the two different 
forms of these solutions are convertible into each other. This paper also is an application 
of several complex variable analysis and its usefulness. In future investigations, one may 
seek the solution of the inverse problem for transport equation with attenuation using this 
several complex variable analysis. 
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